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Abstract 

Within the constituent quark model, we examine the extent to which the 
deuteron photo-disintegration at 2-6 GeV can be described by the quark- 
exchange mechanism. With the parameters constrained by the np scattering, 
the calculated differential cross sections disagree with the data in both magni- 
tude and energy-dependence. The results can be improved if we use a smaller 
size parameter for quark wavefunctions. We also find that the on-shell ap- 
proximation used in a previous investigation is not accurate. 

1. Introduction 

It has been well established that nuclear dynamics at low energies can 
be described in terms of hadronic degrees of freedom. Within the context 
of Quantum Chromodynamics(QCD), this hadronic picture is expected 
to break down at sufficiently high energies where the basic degrees of 
freedom must be quarks and gluons. It is therefore interesting and im- 
portant to know in which energy region the transition from hadronic 
picture to quark-gluon picture takes place. This has been the moti- 
vation of the experiments on deuteron photo-disintegration(7d — ► np) 
conducted at SLAC 1,2 around 1990 and at Jefferson Laboratory 3,4 in 
the past few years. More measurements are expected in the near future. 

The limitation of the hadronic picture of deuteron photo-disintegration 
at energies above about 1.5 GeV was first suggested in Ref.5. The data 
at 1-6 GeV accumulated at SLAC and JLab appear 4 to follow the quark 
counting rule 6 in the kinematic region where the momentum transfer of 
the reaction becomes larger than a critical value. On the other hand, 
one expects that the non-perturbative dynamics, such as those due to 
meson-exchange, still play an important role in the transition region. It 



is obviously a highly non-trivial problem to interpret these data. At- 
tempts in this direction have been made in Refs.7-9. 




Fig. 1. The quark-exchange mechanism of -yd — > np reaction. 

In this contribution we would like to report on our recent investiga- 
tion of the deuteron photo-disintegration at 2-6 GeV. We are motivated 
by the quark-exchange model of Frankfurt et al. 9 Based on the quark- 
exchange mechanism illustrated in Fig.l, the authors of Ref.9 succeeded 
in describing the data at 90 degrees both in magnitude and energy- 
dependence. On the other hand, their predictions at other angles failed 
to describe the data unless a phenomenological factor is included. 

Because of the complexities of the quark-exchange amplitude, several 
simplifications were made in the calculation of Ref.9. First the on- 
shell approximation was used to simplify the loop-integration over the 
quark propagator. Then the -yd — > np cross section is cast, with some 
arguments, into a product of a factor depending on a convolution over 
the deuteron wavefunction and a factor depending on an appropriately 
evaluated elastic np cross sections. No explicit calculation in terms of 
quark-gluon degrees of freedom is needed. While the simplicity of their 
final expression of the differential cross section is very appealing, it is 
possible that the difficulties they encountered as well as the successes 
they achieved could be due to the use of these simplifications. The 
purpose of this work is to address this question within a model in which 
calculations of the quark-exchange amplitude can be performed exactly. 

2. Quark-Exchange Amplitude 

We first observe that apart from the photon-quark interaction vertex, 



the exchange mechanism illustrated in Fig.l is the same as that studied 10 ' 11 
in the constituent quark model of nucleon-nucleon interaction. We there- 
fore will perform the calculation within the same theoretical framework. 
Thus our approach is not fully rclativistic, like all of the previous studies 
of multi-quark systems. But this is the only framework within which an 
exact calculation of the quark-exchange amplitude can be performed. It 
is obvious that our intention here is not a realistic confrontation with 
the data. Rather it is aimed at addressing the theoretical questions 
mentioned above. 

In the center of mass frame, the yd —> np amplitude can be written 
as (suppressing the spin-isospin indices) 

t{p,k)=< x { ~ ] \X(k)\$ d >, (1) 

where k is the incident photon momentum, p the final proton momen- 
tum, $ d the dcutcron wavefunction, and Xp~^ the final neutron-proton 

scattering wavefunction. The operator X(k) represents the quark-exchange 
mechanism. The matrix element Eq.(l) can be written as 

t(p,k) = J dPfX^iPFMP^k), (2) 

with 

T(P F , k) = J dPi J dP 2 S(P 1 - P F ) j dP[ j dP! 2 8{P[ +P! 1 + k) 

X{P u P 2 ,P[,P^k)^ d {^^). (3) 

where P\ and p2{P{ and P 2 ) denote the momenta of the initial(final) 
two thrcc-quark clusters. By assuming that the wavefunctions of all 
three-quark clusters in Fig.l are of Gaussian form (<fr(r) ~ exp{^-) 
in coordinate-space), we can integrate out most of the quark momenta 
analytically and obtain 
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dqi2dq 5( > e 2 e 2 



J dqdqp V{q)G(q p )h{q p ,k) 

exp[-b 2 {a kk k 2 + a kq k ■ q + a kqp k ■ q p + a kPl k ■ Pi 
+a kPp k ■ P F + a qq q 2 + a qqp q ■ q p + a qPl q ■ Pi 
+a qPF q ■ P F + a qpqp q 2 + a qpPl q p ■ Pi 

+ a q P P F <lp ■ p F + apjPjPf + a Pl p F Pi ■ P F 
+a PF p F P 2 F }], (A) 



where a kk = 19/24, a kq = -1/2, a qq = 3/2, a kqp = -5/2, a 9:9p = 
a q P q P ~ 3, a kPi = 3/2, a qP[ = -1, a 9p p 7 = -3, ap T p T = 7/6, a kPp = -1, 
a g p F = 2, a qp ^p I = 3, ap J p F = —2, and ap F p F = 7/6. The parameter 6 
in Eq.(4) defines the range of Gaussian function; 

In Eq.(4), V(q) is a gluon-exchange interaction, G(q p ) is a quark 
propagator, and h{q pi k) — eqU^k^j^u- describes the photon-quark 

coupling. For simplicity, we follow a quark-model study of ttN scattering 14 
to set 

V(q) = £[-™ s (A* • Aj ) (<Ti • ^i], (5) 

where Aj is the color SU(3) generator, a s = 1 and A s = 1.087 GeV were 
determined from fitting the ttN scattering data. The quark propagator 
is taken to be 



G(q P ) = 



1 



W-y/M 6 + <$-y/m? g + q*+ie 
P 



W-y/M 5 + qj-y/m? q +qj 

-m6(W - ^M 5 + q%- yj m \ + ql) (6) 

where P denotes taking the principal-value integration, m q = 330 MeV, 
and M 5 = 5m q . 

To carry out the multi-dimensional integration in Eq.(4), we use 
the method of Ref. 15 to expand the gluon-exchange interaction and 
deuteron wavefunction in terms of Gaussians : V(q) = y~]f v Vie~ aiq 
and = J2^ G 9j e ^ f3jq2 ■ We find that very precise expansions can 

be obtained with Ny = 35 and No = 50. Substituting these expansions 
into Eqs.(3), and (4), we then obtain 

rrm * 1 1 51233 (a- a- 

N v N G , 
i 3 

exp(B) J dqdP/ cxp +s-^j (7) 

where x = (q, Pi) is a two-component vector and 

B = -b 2 {a kk k 2 + a kqp k ■ q p + a kPp k ■ P F + a qpqp qp + a qp p F q p ■ P F + a Pp 

A = 



2b 2 (a qq + p-) a qPl 

a qPl 2b 2 (a PlPl + §■) 

b 2 (+a kq k + a qqp qp + a q p F P F ,a kPl k + a qpPl q p + a Pl p F P F ). 



All integrations in Eq.(7) except that over the propagating quark 
momentum q p can be carried out analytically by using the following 
formula 15 for a multidimensional gaussian 

where d is the dimension of the space and iV is the number of variables, 
x = {x[,x~2, ...,x~n} is a N-component vector. Schematically, Eq.(7) is 
then reduced into 

N v N G 

T(P F , k) = [J2 V "9i / dq P G(q p )h(q p , k)I id (P F , k, q p )\ (9) 

where /^ (Pp, k, q p ) are from integrating analytically various Gaussian 
functions. We then carry out the remaining integration over the propa- 
gating quark momentum q p in Eq.(9) numerically by using the standard 
substraction method to handle the singularity of G(q v ) of Eq.(6). 




. 2. Some results from our fits to np elastic scattering data. The data are from Refs. 12-13. 



3. Results 

The first step to evaluate the jd — ► np amplitude Eq.(2) is to con- 
struct a np model for generating the scattering wavefunction xp It 
can be written in terms of np scattering t-matrix 



where E = 2Epj(p) is the total energy. We have found that the np 
elastic scattering data 12 ' 13 in the relevant 4-12 GeV region can be fitted 
by parameterizing the amplitude t(p',p,E) as a sum of the exchanges 
of it, <t, p, u, and Pomeron. Our fits at four energies are displayed in 
Fig. 2. The details will be given in Ref.16. 
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Fig. 3. The 7ci — > np differential cross sections. No np final state interaction is included. The 
solid curves are from an exact calculation of Eq.(7), while the dashed curves are from taking the 
on-shcll approximation used in Ref.9. The data are from Refs.f-4. 



To explore the dynamical content of the quark-exchange amplitude 
T(Pp, k), we first perform calculations with the final np scattering ne- 
glected; i.e. setting Xp~^*(j7) = $(p ~ P) m evaluating Eq.(2). The 
dcuteron wavefunction is generated from the Paris potential. With the 
specifications given in section 2, the only free parameter of the calcula- 
tions is then the range parameter b of the quark wavefunctions. We find 
that the predicted cross section is very sensitive to this parameter. If we 
set b = 0.52 fm used in the quark model of nuclcon-nucleon interaction 11 , 
the predicted cross sections are an order of magnitude smaller than the 
data. This is the first major difficulty of the quark-exchange model con- 
sidered here. To get reasonable magnitudes at 90° degrees, we need to 
use a much smaller value b = 0.25 fm. The calculated differential cross 



sections are the solid curves shown in Fig. 3. 

We now turn to investigating the validity of the on-shell approxi- 
mation used in Ref.9. This amounts to performing our calculations by 
neglecting the principal- value part of the quark propagator Eq.(6). The 
calculated results are the dashed curves in Fig. 3. The difference between 
the solid and dashed curves indicate very sizable contributions from the 
principal- value parts of the integrations. Our results strongly suggest 
that the on-shell approximation used in Ref.9 is not accurate. 
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Fig. 4. The ■yd — > np differential cross sections. The solid curves are from our full calculations, 
while the dashed curves are obtained when the final np scattering is neglected. The data are from 
Refs.1-4. 



Our complete calculations are the solid curves in Fig. 4. The dashed 
curves arc obtained when the np final state interactions are omitted. We 
see that the np final state interaction, constrained by the data shown in 
Fig. 2, do not change our predictions significantly except at 70 degrees. 
Clearly our complete calculations are not able to describe the data. 

To end, we emphasize that this is a very first step toward a realistic 
calculation of jd — > np reactions in terms of quark-gluon degrees of free- 
dom. To firmly assess the quark-exchange model, several improvements 
must be made in our calculations. First, we find that our predictions 
are sensitive to the choice of the quark mass. This suggests that the use 



of the quark propagator Eq.(6) with a constituent quark mass 330 MeV 
may not be correct. We may need to include the momentum-dependence 
of quark mass, as suggested by some QCD model calculations 17 . Sec- 
ond, the gluon-exchange interaction Eq.(5), taken from a ttN study, 
needs to be improved to include more realistic momentum-dependence 
and other components of the gluon-exchange interaction. Third, we find 
that the main contribution to the cross sections is from the high mo- 
mentum component of the deuteron wavefunction in the region where 
the employed nonrelativistic Paris potential may not be valid. It is nec- 
essary to generate the deuteron wavefunction from a model which can 
also describe np scattering data at high energies, such as those shown 
in Fig. 2. It may be possible to obtain such a model by extending the 
ttNN model of Ref. 18 to higher energies. Finally, the size parameter b 
of the three-quark clusters needs to be chosen from a more fundamental 
consideration. 
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